Fragmentation pathways of nanofractal structures on a surface 
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We present a theoretical analysis of the post-growth processes occurring in 
nanofractals grown on a surface. For this study we have developed a method which 
accounts for the internal dynamics of particles in a fractal. We demonstrate that 
the detachment of particles from the fractal and their diffusion within the fractal 
and over the surface determines the shape of the islands remaining on a surface af- 
ter the fractal fragmentation. We consider different scenarios of fractal post-growth 
relaxation and analyze the time evolution of the island's morphology. The results 
of our calculations are compared with available experimental observations, and ex- 
periments in which the post-growth relaxation of deposited nanostructures can be 
tested are suggested. 
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I. INTRODUCTION 



Nowadays, nanoscience is a rapidly developing research domain This generic word 

refers to the study performed on systems having a characteristic length scale on the order of 
a nanometer: a length scale at which new specific physical and chemical properties emerge 
in the system. One of the main goals of nanotechnology is the development of controlled, 
reproducible, and industrially transposable nanostructured materials jl-7|. In this context, 
controlling of the final architecture of such materials by tuneable parameters is one of the 
fundamental problems. 

The conventional technique of thin- film growth by deposition of atoms 2|, |3| , small atomic 
clusters js, 4| and molecules [l, 2, 8-lo| on surfaces gives a possibility to construct materials 
with pre-defined properties. Recent experiments show that patterns with different mor- 
phology can be formed in the course of cluster deposition on a surface jsl, 0, [ll|. Among 
other possible shapes, droplet-like and fractal islands have been observed in various systems 



otner DO 



13| . It was shown that the island morphology depends on various factors, such as 



the temperature 
roughness [19 



l4, list], particle size 16|, particle deposition rate Is, [iTI, [l8|, substrate 



14 




2l| and interparticle 



20| . concentration of impurities in the system 
interaction energies 0, S]. It was also demonstrated that the patterns on a surface strongly 
depend on the type of the substrate. For example, experimental studies of silver clusters 



depositee on si icon at room temperature showed that droplet-like islands are formed 22|, 



while in []J, [15 



21| it was demonstrated that dendritic shapes emerge on graphite. 



The investigation of the dendritic structures (fractals) has attracted considerable atten- 
23-27|. The formation of such systems provides a natural framework 



tion (111, M 



15 



21 



for studying disordered structures on a surface because fractals are generally observed in 



ar-from-equi 



Fe-N |26|] clusters and Ceo molecules [29 



ibrium growth regime. For example, fractals consisting of Ag |l4 . 

3fl 



15 



2l|, Au 



with the use of the cluster deposition technique 



lave been fabricated on different surfaces 



29 



The growth process of fractals has been extensively studied in experiments 17H20 



3l|. In 



|l7,[l 



28|, 



L8[ a quantitative experimental study of spherical antimony cluster diffusion 



on graphite was performed. It was shown that the size of the emerging fractals depends on 
the cluster deposition rate. The influence of cluster size on fractal morphology was experi- 



mentally studied in 



16| . In that work antimony clusters of different size were subsequently 
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deposited on graphite surface, and it was demonstrated that the fractal branch width de- 
pends on the size of the deposited clusters. Molecular processes, underlying the Cgo-fractal 
formation on graphite substrate were investigated experimentally by use of the scanning 
tunneling microscopy {29 1. The self-organization of silver clusters on graphite surfaces with 
different crystallographic orientations was experimentally investigated in loj]. It was shown 
that the size of the formed fractals depends on the crystallographic planes of graphite, which 
influences the cluster mobility over a surface. 

Contrary to the process of fractal formation, the process of post-growth relaxation and the 
question of stability of deposited structures are still not well understood. The understand- 
ing of the post-growth relaxation processes would allow one to control the self-organization 
processes of particles on a surface for the purpose of obtaining patterns with predictable 



morphology. An illustrative example of pattern manipulation was given in 31| by adding 



metal impurity to the system. In that work different morphologies of Cqo films with triangu- 
lar, dendritic and fractal-like (lll)-oriented single-crystal grains were detected by changing 
the thickness of the pristine fullerene film and the concentration of Ag impurities. 



The post-growth transformation of si 
surface was experimentally studied in 



ver cluster fractals to compact droplets on graphite 
[isi 21 1 . It was demonstrated that depending 



on the experimental conditions the shape and the size of the stable silver droplets changes 
significantly 21 1. In jl^. Iisl. 21 1 it was shown that oxidizing of silver clusters results in rapid 
fragmentation of a fractal, leading to the formation of several compact droplets. 

An important characteristic, which determines fractal formation and post-growth relax- 
ation dynamics is the mobility of a cluster on the substrate, which in turn is temperature 
dependent. Fractals of gold clusters, grown at room temperature on ruthenium substrate 



32| . Thermal 



fl^ . I21I . In these papers 



undergo a transformation into compact droplets after annealing at 650 K 
relaxation of silver cluster fractals was experimentally studied in 
it was demonstrated that due to thermal annealing the fractal branch width increases and 
eventually the fractal breaks into smaller parts. 

The dynamics of particles on a surface was also studied theoretically. An efficient theo- 
retical tool for describing particle dynamics on a surface is the diffusion limited aggregation 
(DLA) method 33|]. In this method each particle on a surface moves freely in a random 
direction until it collides with another particle. In the case of collision both particles stick 
together and become immobile. The DLA model was used for a qualitative description of 



34|. 



the process of fractal formation on a surface 

Two-dimensional theoretical model based on the DLA method has also been developed 

n 

for the description of thermal relaxation of fractals on surfaces [3^ . In this model particles 
are treated as immobile only if they are complectly surrounded by other particles. In all 
other cases particles are allowed to move along the branches of a fractal with a certain 
probability |3J|. The parameter dependent method developed in [3J] was used to describe 
the thermal transformation of a fractal into a droplet. The description of fractal instability 
in 3^ was limited only to one particular choice of parameters. However, the correspondence 
of these parameters to the actual experimental values was not established. 

The island size distribution function is a fundamental quantity in the kinetic description 
of island growth. It has been widely used to characterize the experimentally measured 
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37| as well as computed (38H42[ surface morphologies. The scaling of the island size 



distribution, island density, monomer density and other morphology characteristics of the 



system are often studied function of the sur 



ratio of the diffusion rate to the deposition rate 



ace coverage, and as a function of the 



38|-|4l| . The scaling of the afore mentioned 



quantities allows to determine important physical parameters describing the kinetic growth 



processes on a surface, as, e.g., the activation energy and the diffusion constant 
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4ll- 



An important characteristic of the island size distribution scaling is the scaling function. 



The shape of the scaling function is independent of the initial 
surface and is determined by the island's growth mechanism 



distribution of particles on a 



43| . The scaling of the island 



size distribution emerging during the post-deposition growth processes on a surface was 



also suggested 



44| . The scaling of island morphology characteristics has been performed 



for various model systems 38H41|, and allows to characterize complex kinetic processes 
occurring on a surface. 

In the present paper we make an important step towards understanding of the differ- 
ent evolution scenarios of nanofractals morphology. We present a systematical theoretical 
analysis of the post-growth processes occurring in nanofractals on a surface. For this study 
we have developed a method describing the internal dynamics of particles in a fractal with 
accounting for their diffusion and detachment. We demonstrate that these kinetic processes 
determine the final shape of the islands on a surface after the post-growth relaxation. We 
consider different scenarios of fractal relaxation and analyze the time evolution of the island's 
morphology. The results of our calculations are compared with experimental measurements 



of the post-growth relaxation of silver cluster fractals on the graphite substrate [isl, |21 |. 
In particular, we analyze the island size distributions calculated at various conditions and 
different post-growth fragmentation regimes. In conclusion we outline a number of open 
problems which should be investigated in the future. For instance as a possible next step 
one could explore the scaling of the island size distributions during the fractal fragmentation 



processes, as also suggested in 4J]. This and many other interesting relevant questions are 



beyond the scope of the present work and are left for further investigations. 



II. THEORETICAL METHODS 

In this section we discuss the theoretical methods used for studying the dynamics of 
particles on a surface. Computations were performed with the use of the MBN Explorer 
computer package, which is developed for structure optimization, simulation of dynamics 



and growth processes in various nanosystems 45l-l55l|. Below we describe the general idea of 



the computational method used in our work and explain how internal dynamics of particles 
in a fractal has been accounted for. 

To study the diffusion of particles over a surface we used a version of the kinetic Monte 
Carlo (KMC) method [3, IsgUs^. The KMC method is based on the Monte Carlo algorithm 
and is widely used for the study of time evolution of various processes occurring in nature 



57H63|. such as surface diffusion of p articles |64| . vacancy diffusion in alloys j65|, damage 
accumulation and amorphization j66| and many others. The processes described using the 
KMC method always occur with certain predefined rates. Note, that these rates are input 
in the KMC algorithm, and the method itself cannot predict them. The calculation of the 
kinetic rates for different processes is usually a nontrivial problem. The kinetic rates are 
material-dependent parameters of the KMC method, which in the case of particle diffusion 
over a surface are determined by the atomic composition of particles, substrate material 
and interparticle interactions. Therefore, by varying values of the kinetic rates, the KMC 
method can be used to study the dynamical behavior of various molecular systems. Advanced 
computational methods are often necessary for the kinetic rates calculation. For example, the 
rate of particle diffusion over a surface can be extracted from molecular dynamics simulation 



(for a review see, e.g., |67|). 



The idea of the KMC method is as follows. The time-evolution of a molecular system 
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is modeled stepwise in time. With a certain probability, at each step of the simulation, 
the system undergoes a structural transformation. The new configuration of the system is 
then used as the starting point for the next simulation step. The transformation of the 
system is governed by the kinetic rates, input into the KMC method. Note, that at each 
simulation step the system can be transformed into one of several states. Thus, in the KMC 
method, the probability for the system to attain a certain configuration is proportional to 
the corresponding kinetic rate. Due to its probabilistic nature, the KMC method allows 
to study dynamical processes on time scales significantly exceeding the time scales of the 
conventional molecular dynamics simulation. This method is ideal in the situations when 
the intermediate details of the dynamical processes are not so essential, and the transition 
to the final state of the system can be parameterized by a few kinetic rates. 

In our studies we used a modification of the conventional KMC method. First, for each 
particle on a surface we determine the number of possible diffusion directions. Thereby, a 
particle can either diffuse freely over a surface, or diffuse along the periphery of the already 
pre-formed structure on a surface. A diffusion direction of the particle is chosen randomly 
in such a manner that all possible diffusion directions are equally probable. For a given 
diffusion direction the probability of particle diffusion is then calculated, and the particle 
is moved in this direction according to the calculated probability. Below we discuss this 
method in detail. 



A. Fractal Growth 



To model the growth of a fractal on a surface we used the diffusion limited aggregation 
(DLA) method 33|]. Using a module of the MBN Explorer program 68| we have computed 



the growth process of a fractal by depositing particles on a surface in the vicinity of the pre- 
defined growth center. To compare with the experimental measurements 2l| we have 



used in our simulations the model parameters consistent with the experiment. Thus, the 
diameter of a particle has been taken to be 2.5 nm, which corresponds to the size of an Agsoo 
cluster used in Q, [l5, 21|. The deposition fiux has been decreased from Fstart = 7.2 x 10^^ 
particle/cm^s to Fend = 1-1 x 10^^ particle/cm^s because the area to which the particles 



are added decreases as the size of the fractal increases. The used values of the particle 



deposition fiux are chosen higher than the experimental value reported in [l^ [2l|, F 



FIG. 1. Island of deposited particles on hexagonal grid. Processes essential for a fractal formation 
on a surface are shown by arrows, T is the free particle diffusion rate, Td{m, n) is the diffusion rate 
along the periphery of fractal or island, and re(/) is the particle detachment rate from fractal or 
island. The diffusion rate along the periphery depends on the number of broken bonds (m) and 
the number of maintained neighboring bonds (n). The particle detachment rate depends on the 
number of broken bonds (/). In the depicted example m = 3, n = 1, I = 2. 

particle/cm^s, in order to accelerate simulation of the fractal growth. The simulated fractals 
have been used as initial structure in the investigation of fractal fragmentation. To speed up 
the calculation, we simulated particle dynamics on a 2D hexagonal grid, on which a particle 
has up to six neighbors, as illustrated in Fig. [H The size of a single grid cell in this case is 
defined by the particle diameter. 

To simulate fractal growth the following procedure has been adopted. At every step of 
the simulation new particles are deposited on the surface according to the deposition rate 
and occupy some of the free cells in the grid. Simultaneously, the already deposited particles 
diffuse over the surface with the rate 
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kT 



(1) 



r = z/i exp 

where Ea is the activation energy, ui is the attempt escape rate, T is the temperature of the 
system and k is the Boltzmann constant. The process of particle diffusion over the surface 
is schematically illustrated in Fig. [1] 

An important quantity in the DLA method is the time step. At, which defines the char- 
acteristic time for particle diffusion over a surface as 



At = 1/r. (2) 

The time step At is related to the coefficient D of particle's diffusion over a surface arising 



in the equation of diffusion 



69 



70|. The solution of the diffusion equation in two dimensions 



gives the probability to find a particle at the instant t being at distance [r, r + dr] from its 
initial position as 



u{r, t)dr 



exp 



rdr. 



(3) 



2Dt \ ADt^ 

Jsiiig this probability function one derives the mean square displacement of a particle as 



69 



7Q| 



POO 

(n - ro)2 = / u;(r, ti - t^ydr = 4D{t, - to), (4) 

^ 

where tq and ri are the distances to a particle from the initial position at two successive 
instances to and ti. Equation (jlj) allows to express the diffusion coefficient as 



D 



(Ar^) 
2zAt ' 



(5) 



where (Ar^) is the mean-square displacement of a particle per time At, and z is defined by 



69 



70|. In the case of particle diffusion over a surface z = 2 



the dimensionality of space 
(see Eq. (HD). 

On the other hand the mean-square displacement depends on the diffusion rate and on 
the particle hopping length, which is in the example considered is equal to the particle 
diameter do: 



(Ar^) = TdlAt. 



(6) 
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Here At has a meaning of a single simulation step defined in Eq. ([2]). Substituting Eq. ([6]) 
into Eq. ([5]), one derives 

D - ^. (7) 
Equation (I7j) allows to estimate F (and therefore At) once the diffusion coefficient is known: 

At = (8) 
2zD ^ ' 

The diffusion coefficient of an Agsoo cluster on graphite at room temperature was mea- 
sured as 2 ■ 10"'^ cm^/s [2]J. Substituting this value into Eq. (jHl), one obtains At = 78 ns. 

Substituting Eq. ([1]) into Eq. ([7]), one relates the diffusion coefficient to the activation 
energy and temperature: 

It follows from Eq. (Q that the diffusion coefficient decreases as the activation energy grows. 
This results in an exponential growth of the time step At with Ea, since At ~ 1/D (see 
Eq. (IHl)). Equation ([8]) introduces the optimal time step for the computations, because it 
defines the characteristic time at which a freely deposited particle gets displaced for the 
distance do, i.e. to the neighboring lattice cell (see Fig. [1]). Note, that the diffusion of 
deposited particles is typically the fastest process in the system. Therefore, it determines 
the minimum time step and the time scale for the growth and fragmentation of fractal 
structures on a surface. It is computationally inefficient to perform simulation with time 
steps much less than At, because in this case the particles will practically not move during 
t < At. 

In our method we have employed the following procedure to model particle dynamics over 
a surface: for each step of the simulation the deposited particles without neighbors have six 
possibilities for diffusion, see Fig. [TJ The direction of displacement is chosen randomly and 
each particle is moved to a neighboring lattice cell in the chosen direction. Thereby each 
step of the simulation corresponds to At = 78 ns, as estimated above. The particles at the 
fractal periphery diffuse slower, as shown in the next section, and therefore are displaced 
less frequently. 
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B. Kinetic processes in fractal fragmentation 

In this work we consider fragmentation of a two-dimensional fractal consisting of rigid 
particles of the equal radius. The relaxation of such fractal is driven by the diffusion of 
particles along the fractal periphery and by the detachment of particles from the fractal. 
Both processes are schematically depicted in Fig. [H The diffusion and the detachment 
rates depend on the activation energy and the particle-particle interaction. In the Arrhenius 
approximation the diffusion rate of a particle along the fractal periphery reads as: 



Tdim, n) = z/2 exp 



mEh nAe 



kT 



(10) 



kT kT 

where m is the number of bonds that are broken due to the particle motion, Eb > is the 
binding energy between two particles, n is the number of maintained neighboring bonds 



between two particles and Ae < Ei, is the diffusion energy barrier [15|, |3J| , 1^2 is the attempt 
escape rate. Equation (HU]) describes the probability of a particle to overcome a potential 
energy barrier, which for a particle diffusing along the fractal periphery is parameterized by 
the energies Ei,, Ae, and Ea- Note that the parameter Ea, which enters Eq. f lTOj) . depends 
on the simulation time step At, as discussed in the previous section. Therefore, only the 
parameters Eh and Ae define the potential energy barrier for particle diffusion along the 
fractal periphery, while Ea characterizes the time scale. 

Note that Eq. f lTU]) does not account for the bonds which may be created in the system 
when a particle diffuses. This feature of Eq. ffTOj) is easy to understand. Indeed, the particle 
diffusion process in our model is considered stepwise, i.e., at each step of the computation 
a particle is displaced with a certain probability in a random direction for the distance 
equal to its diameter. But prior the particle is displaced to its new position there is no 
information about the newly created bonds in the system (causality principle). Therefore 
only those bonds which the particle forms with its neighbors prior the displacement influence 
the diffusion dynamics in the system. 

The evaporation (detachment) rate of a particle from the fractal is given by 



'kT 



A/i 



kT 



where / is the number of bonds broken during the particle detachment 



chemical potential change associated with particle detachment 



15 



Tom fracta 

g 



34 



(11) 

A/i is the 



72| . is the 
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attempt escape rate of a particle in its equilibrium state. Equation (fTTj) can be understood 



71| . which studies the liquid-H-gas 



within the framework of the classical nucleation theory 
transition in droplets. It is written in the Arrhenius approximation, similarly to Eq. (ITOl) . 
For the further discussion of the fractal fragmentation we put: 



1/2 ~ 1/3 = T^- (12) 

Such situation occurs when the characteristic attempt escape rate of a particle leading to 
its diffusion or detachment are close. Equations ([TOl)-( |TT|) describe the dependence of the 
probability of different essential kinetic processes on the values of Ea-, Ef,, Ae, A/i, which 
below are called the model parameters. For convenience, in this paper we define all the 
model parameters in units of kT (1 kT = 0.026 eV) at room temperature (300 K). 



III. RESULTS 



A. Model parameters 

The interaction energy between the deposited particles and the substrate is responsible 
for the particle mobility over a surface, as follows from Eq. ([1]). The energetic parame- 
ters for the atomic-scale processes on the Ag(lOO) surface were studied in [73|, while the 
energetic parameters for Si atom migration on Si(100)-2xl surface were discussed in 56 |. 
The interaction energy of Agsoo (E^^), Ceo and Sb2300 5 (E^^) clusters with graphite 



7^, = 6.9 kT [75| 



surface at room temperature has been estimated as E^^ = 6.6 kT 

and E^^ = 27.1 kT 17[. The significant spread of the values indicates the essential role of 
interatomic interactions in defining the activation energy. The value of Ea defines the time 
scale of the fractal growth and fragmentation processes, as discussed in section III A[ In this 
paper we describe the dynamics of silver cluster fractals with the above mentioned value of 
Ea- Another important quantity characterizing the particle diffusion over a surface is its 
attempt escape rate u (see Eqs. f[T0|) - f|T2|) ). which can be estimated as 



2Dz 



exp 



E., 



kT 



(13) 



For a silver nanoparticle with do = 2.5 nm deposited on graphite the diffusion coefficient at 
room temperature D ~ 2 ■ 10^'' cm^s~^ 21|, resulting in u 



0.94 ■ 10^° s~^ Note that the 
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attempt escape rate for silver clusters is 2-3 orders of magnitude smaller than for individual 
atoms on a surface [3, S, 0, Q • 

The interaction energy between two particles, Eb, depends on the atomic composition of 



I21I I . Thus, it was shown 



the particles and on the presence of impurities in the system [lJ.[l5. 
that the presence of oxygen impurities in a silver cluster deposited on graphite leads to the 
decrease of and consequently to the reduction of fractal stability. A systematic study 
of the afore mentioned factors on the interparticle interaction energy is beyond the scope 
of this paper and deserves a separate systematic investigation. Here for us is important 



that according to experiment [IJ, l2l| silver cluster fractals are formed and decay on the 
comparable time scales. This means that in our model Eb should be of the same order 
of magnitude as Ea- It is worth noting that the time of fractal formation can become 
significantly smaller than the time of fractal fragmentation if the conditions at which the 
fractal is kept after growth rapidly change, e.g. the temperature of the system is increased. 
The diffusion barrier energy Ae depends on the atomic com pos ition of the cluster and usually 



76| 



amounts 0.05 — 0.2 of the bonding energy of two clusters p 

The change in the chemical potential Afi arises due to the energy difference caused by the 
change of the number of particles in the system. The chemical potential characterizes the 
ability of particles to diffuse from regions of high chemical potential to those of low chemical 



potential and is defined as the partial derivative 



where U and S are the total energy and the entropy of the system, V is its volume and 
is the number of particles in the system. The variation of the chemical potential arising 
due to a structural transformation in the system can be calculated from the known values 
of the chemical potential of individual components of the system before and after the trans- 
formation. For example, for the evaporation of a silver nanoparticle from a fractal with 
particles on graphite surface 



Ag(fractalN) + C(graphite) — )■ Ag(fractalN-i) + C(graphite) + Ag(particle), (15) 

the corresponding change of the chemical potential can be calculated as a difference between 
the chemical potential of the products and the educts. With AiAg(fractaiN) ~ /^Ag(fractaiN-i) 
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obtains 



A/i — /iAg(particle) • 



(16) 



The chemical potential 



can be measured experimentally 78| and is tabulated for many 



substances (see e.g. [79|,l80[). It depends on the phase state of the system: for the gas of silver 



79|. These 



atoms /iAg'^'' = 2.55 eV, while for the silver in the liquid phase yuj^'g^"^'' = 0.8 eV 
values and Eq. (fT6l) allow one to suggest that the change of the chemical potential in the 
silver fractal fragmentation process, at room temperature lies within the range 30 — 100 kT. 



B. Fractal growth 



We stress that this paper is devoted to the process of fractal post-growth evolution. 



The problem of nanofractal formation has been studied in many papers (see e 



20 



28|, 



29 



31 



17 



34j ) and we intend to extend this analysis in a separate publication. Therefore, 



here we give only a hint how the fractals of the morphology of interest can be obtained. 




4- 



^ 200 nm 




500 nm 



FIG. 2. (a) Fractal structure simulated using the DLA method; (b) structure of silver cluster 
fractal experimentally grown by clusters deposition technique on graphite surface [l^; 



Using the method described in Sec. IIA we have simulated several fractals with the struc- 
ture being very similar to the silver cluster fractals grown experimentally on graphite surface 
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Q,|l5|, 



2lj . The fractal structure shown in Fig. [2^ have been chosen for the further inves- 
tigation of the post-growth relaxation processes in fractals. The diameter of the fractal is 
635 nm, which is somewhat smaller than the diameter of the experimentally grown struc- 
211 . For the sake of illustration in Fig. [2b we show the experimen tally gr own 



tures 
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15 



silver cluster fractal prior thermal annealing, which triggers fractal fragmentation 14l.ll5l.l21| 



The important characteristic of a fractal 



fractal dimension is generally defined as 81 



is the fractal dimension df. The Hausdorff 



82|: 



(17) 

1^0 log[l//] ^ ^ 

Here N{1) is the number of self-similar structures of linear size / needed to cover the whole 
structure. In practice the fractal dimension is usually calculated by the box-counting method 
Equation ( 1T7|) has been used to calculate the fractal dimension of the structure shown 



in Fig. [2^. This calculation resulted in d^j^ = 1.76. This value is in a good agrement with 
experiment for silver cluster fractals grown on the graphite surface, which gives d'^^^^ = 
1.7±0.1 Q. 

As illustrated in Fig. [2] the topology of fractals simulated by the DLA method is very 
similar to the fractal topology grown in experiment. In both cases the fractals shown in Fig. [2] 
have several main branches, growing from the center of the fractal. The branch width of the 
fractal simulated by the DLA method is ~ 10 nm, while the typical experimental width of 
the branch is 15 — 30 nm [l^. The difference arises because the particles in the simulation 
were deposited on a surface at a higher rate than those in experiment (see section III Al for 
details). In addition, the implemented DLA method does not allow the deposited particles 
to be placed atop of a growing fractal. Another factor affecting this difference is the sticking 
probability of the deposited particles assumed to be equal to one, meaning that if a particle 
meets another particle on a surface the two particles stick and do not move together. This 
is probably not the case in experiment, where the sticking probability can be lower than one 
and the mobility of complexes with two or more particles is no equal to zero. The sticking 
probability is less important in the fractal fragmentation process as this process is mainly 
governed by the detachment rate Fe, and the diffusion rate F^ introduced in section III Bl 
Since the main focus of this paper is the investigation of pathways of fractal fragmentation, 
we do not discuss further the effect of sticking probability on pattern formation. 




FIG. 3. Evolution of fractal structure on a 650 x 750 nm^ substrate with periodic boundary condi- 
tions. The initial fractal structure shown in the middle undergoes fragmentation in different final 
states depending on the interparticle interactions in the system. Numbers above the corresponding 
images indicate the values of Ef, and Ae used in the simulations (in units of kT), Afj, = 2 kT in all 
cases. The simulation time t is given for each path of the fragmentation. 

C. Fractals fragmentation 

In this section we perform analysis of the fractal post-growth relaxation using the method 
described in Sec. IIIBi According to the estimates performed in Sec. Ill Al a single time step 
in our calculation is equal to At = 78 ns, which allows one to simulate the process during 
the time period 
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t = Nstep^t, (18) 

where Ngtep is the number of simulation steps. 

In the present work we analyze several paths of fractal fragmentation. The rate of fractal 
decay depends on the interparticle interaction, and it defines the morphology of the frag- 
ments that are formed during the process. Snapshots of the structures arising at different 
stages of the fragmentation process simulated at different parameters of interparticle inter- 
actions are shown in Fig. |3l This example shows how different can be the fragmentation 
paths and the fragments morphology. 

Figure |3] shows that for Ef^ = 1 kT, Ae = 0.2 kT one observes an entire defragment at ion 
of a fractal, which is the fastest fragmentation path. In this case the interaction energy 
between the particles is relatively weak and the probability to evaporate a particle from the 
fractal is much higher than the probability of newly deposited particles to nucleate. This 
fragmentation scenario can be realized in experiment if the temperature of the system is 
rapidly elevated after the fractal was created. 

Figure [3] shows that for > 2 kT the fractal melts in a number of compact droplets. 
Depending on the energies of interparticle interactions the shape of the droplets becomes 
different. Thus, for Eb = 2 kT, Ae = 0.4 kT three large, almost spherical, droplets of a sim- 
ilar size are formed. In this case the binding energy Ef, between the particles is rather small, 
allowing relatively easy detachment of particles, but at the same time it is large enough to 
make the characteristic particle detachment time comparable with the characteristic parti- 
cle nucleation time, thereby preventing the system from entire defragment at ion, observed at 
Eb = 1 kT. Thus, the fragmentation path at Eh = 2 kT goes via the rearrangement of the 
entire system, and the formation of large stable droplets. 

A further increase of the interparticle interaction energy leads to the change of the fractal 
fragmentation pattern. As seen in Fig. [3] at = 4 — 6 kT the fractal fragments into 
several compact droplets. The analysis of morphology of the created patterns leads us to 
the following main conclusions: (i) the growth of Eb leads to the increase of the number of 
droplets on a surface (see Eb = 4: kT and Eb = Q kT) and to the decrease of their average size. 
This happens because the detachment of particles from the fractal becomes energetically an 
unfavorable process, and the fractal fragments mainly due to the peripheral diffusion of 
particles, initiated at the peripheral defect sites, (ii) The increase of the peripheral diffusion 
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barrier energy Ae suppresses the diffusion of particles, resulting in a slower evolution and 
fragmentation of the fractal shape. It is remarkable that at Ef, = Q kT and Ae = 1.2 kT 
one observes the formation of elongated islands on a surface which follow the direction 
of the fractal branches. A further increase of the interparticle binding energy with the 
simultaneous lowering the barrier energy for the particle peripheral diffusion favors the 
formation of elongated islands on a surface. Figure |3] illustrates this for Ef, = 12 kT and Ae = 
1 kT. In this case the timescale for the particles to detach from the fractal is significantly 
larger than that for the peripheral particle diffusion. 

A simultaneous increase of the interparticle binding energy and the barrier energy for the 
particle peripheral diffusion leads to the growth of the fractal life time. Figure [3] shows that 
for Eh = 24 kT and Ae = 12 kT the fractal has no noticeable changes in its morphology 
after 4 s of simulation. In the case when the interparticle energies are large, the fractal 
fragmentation is expected to occur on a larger time scale and can be simulated numerically 
if the value of the simulation time step is increased. 

The important characteristic of the fractal fragmentation is the number of fragments 
at a given time. The smallest fragment is a single particle. The time evolution of the 
number of fragments calculated for different sets of model parameters is shown in Fig. 
and Fig. HJi. Curve 1 in Fig. HIi shows the time evolution of the number of fractal fragments 
at Eh = 1 kT. The number of fragments in this case rapidly approaches the asymptotic 
value, approximately equal to the half of the total number of particles in the fractal. This 
means that the system dominantly consists of dimers. With increasing Eh the number of 
fragments at the equilibrium decreases, as seen in Fig. [31 It is interesting to note that 
at Eh = 2 kT there are three dominating large islands (see Fig. [3]). The total number of 
fragments at the end of the simulation in this case is equal to 100, being much smaller than 
the total number of particles in the system. This feature arises in the situation when a large 
number of single particles detach from the large droplets but later stick back. In this case 
the number of single particles fluctuates rapidly resulting in the large fluctuations of Nfr{t) 
dependence shown in Fig. by curve 2. These results have been calculated for a fractal on 
a 650 X 750 nm^ substrate with periodic boundary conditions. 

Figure mi shows that there is no dramatic change in Nfr{t) dependence with the growth 
of Ae at a constant value Eh. This analysis shows also that the growth of Ae preventing 
particles peripheral diffusion hinders the fast transformation of droplets into compact islands 
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simulation time (s) 

FIG. 4. Time evolution of the number of fragments iVj,., {Rmax) introduced in Eq. (|19p and 
of the (S/P) ratio introduced in Eq. ()27p calculated for the fractal structure shown in Fig. [3l 
The fractal fragmentation have been analyzed at Af.t = 2 kT for different values of the binding 
energy and the barrier energy Ae. Plots (a) — (c) show the results of calculation obtained 
for Ae = 0.2E{, and the different values of the binding energies between two particles. Lines 1-6 
correspond to E'f, = (1, 2, 3, 4, 5, 6) kT, respectively. Plots (d) — (/) represent the results obtained 
at Eb = 4 kT for different values Ae = (0, 0.4, 0.8, 1.0, 3.2, 4) kT. The direction of growth of 
Ae is shown in these plots. 

which eventually results in the increase of the number of fragments on a surface. 

As seen in Fig. [3l in the course of fractal fragmentation the mobile particles can coa- 
lescence into islands, i.e. groups of particles bound together. The size and the number of 
islands on the substrate depend on the binding energy Ef, and the barrier energy Ae. The 
important characteristic of the fragmentation pattern on a surface is the average maximal 
radius of the created islands which reads as 

1 

^^fr i=l 

where Nfr is the total number of islands on a surface, R^ax is the maximal radius of the i-th. 
island. The dependencies of {Rmax){t) calculated at different values of and Ae are shown 
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in Fig. HJd and Fig. H^. These figures show that in average (Rmax) approaches the equihbrium 
value at the chosen values of model parameters except for Eh = 3 kT, Ae = 0.6 kT when the 
large fluctuations of (Rmax) develop and grow with time. This happens because at Eb = 3 kT 
the rate of single particle detachment turns out to be so that only several particles are able 
to overcome the detachment energy barrier at one simulation step. The escaped particles 
freely diffusing over the surface after a short period of time return to the same or some other 
island. Although the number of fluctuating fragments on the surface in this case is relatively 
small (see Fig. and Fig. [5]), the fluctuations of {Rmax) become considerable because at 
these conditions small islands can be spontaneously created but most of them disappear just 
after several simulation time steps. Thus, for example, for ti = 3.484 s (i?rn.aa;)i=21.5 nm, 
while for t2 = 3.485 s (-Rmaa;)2=33.4 nm. The change of the maximal radius A{Rmax) in this 
case is 11.9 nm within 1 ms time interval. This happens because for the time frame ti there 
are A^^.^'' = 6 single particles and N^^ = 10 fragments of a larger size with approximately 
equal diameter on the surface. For the time frame t2 the number of large fragments is , 

(2) 

still equal to 10, while there are no single particles on the surface (i.e. Ng = 0). With 
R^^^ = R^^^ = Rl being the characteristic radius of the large island, R^^ = R^P = Rs the 
radius of a sing le particle, and N^^^ = A^f ^ = A^^, one derives 

mmax) = [Rl - Rs) , (20) 

where AA^s' = A^^.^^ — A^^.^'' is the change of the number of single particles. A"! = A'^ + A^^.^^ 

(2) 

is the total number of particles at instance ti, and N2 = + Ng is the total number of 
particles at instance ^2- Substituting values for AA^s, A'l, A"! and N2 in Eq. (!20!) for the 
special case considered one obtains 

A{Rmax) = ^iRL- Rs) ■ (21) 

Substituting {Rl)=32 nm and Rs = 1.25 nm in Eq. fl2T]) . one derives A{Rmax) = H nm. 
Equation (pTj) shows that A{Rmax) increases with R^ which grows with time until it reaches 
the equilibrium value. Equation (|20|) can also be rewritten as 

AA'^A'i 

^^^-"^^ = N! (1 - ANs/N,) - ' ^^^^ 
which shows that for ANs ^ A'l the fluctuation of the average radius A{Rmax) can be 
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Eb (kT) 

FIG. 5. Dependence of Nfr (squares, left scale) and (Rmax) (dots, right scale) on the binding 
energy Ef, calculated for the barrier energy Ae = Q.2Eh after 4 s simulation, corresponding to the 
dependencies shown in Fig. and Fig. [Da. 

several times larger than the the value of the average radius. Note that although the largest 
islands are observed for the model parameter Eb = 2 kT (see Fig. [3]), the largest average 
maximal radius is expected for Eh = 3 kT as depicted in Fig. [51 This happens because the 
number of single particles on the surface for Eb = 3 kT is about 10, while for Eb = 2 kT it 
is exceeding 100. 

Figure ^ shows some dependence of (Rmax) on Ae. The growth of Ae leads to the 
decrease of (Rmax), which is a natural result of a lower peripheral mobility of particles. 

Figure [6] shows the distributions of island sizes in the system after 4 s of simulation. In 
order to improve the statistics, the distributions shown in Fig. [6] have been averaged over a 
time interval r = 0.78 s as follows 
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FIG. 6. Distributions of island sizes formed on the substrate after 4 s of simulation. The distribu- 
tions were calculated at the fixed values of = 4 kT, and A/x = 2 kT for different values of Ae 
as indicated. 



1 fr/2 

Nfr{t) = - J ^ ^.fri^ - (23) 

The histograms in Fig. |6]have been calculated with different barrier energies. The maxima 
in the distributions show the most abundant island sizes. Figure [6] shows that the sizes of 
the islands created in the fractal post-growth fragmentation process depend strongly on the 
binding energy Ei, and the barrier energy Ae. At some values of Ei, and Ae one can identify 
two maxima in the island size distributions. Especially clear this feature manifests itself at 
Ae = kT , Ae = 0.4 kT and Ae = 3.2 kT. The presence of two maxima in the island 
size distributions tells that there are two groups of islands on the surface having different 
preferential island size. 

Let us also analyze the time evolution of the distributions shown in Fig. O Figures [7] 
and E] show distributions of the island sizes calculated at different fragmentation stages 
for a fixed set of the model parameters. Figure [7] illustrates the evolution of the island 
size distribution simulated aX E}^ = 2 kT and Ae = 0.4 kT. After fast fragmentation 
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FIG. 7. Size distributions of islands calculated at different stages of the fractal fragmentation (see 
Fig. S]) for Efj = 2 kT, Ae = 0.4 kT and Afj, = 2 kT. The corresponding simulation time is given 
in the insets to the plots. 



of the fractal into a subset of noncompact islands which occurs on the time scale greater 
than 0.04 s, the distribution of islands sizes has a Gaussian-like shape with the maximum 
centered at 25 nm. In the course of the fractal fragmentation process the magnitude and the 
position of the maximum of the distribution change, because the morphology of the system 
changes due to the evaporation of single particles from the islands and the nucleation of 
single particles. Figure [7] illustrates that small islands nucleate into larger droplets resulting 
in a shift of the maximum of the distribution towards larger island sizes. Interesting that 
the fragmentation/nucleation dynamics leads in this case of study to the formation of two 
maxima which correspond to the presence in the system of the droplets of different radii. 

Figure [7] shows the evolution of the fractal fragmentation process. The initial fragmen- 
tation of the fractal is very rapid. It involves the rearrangement of single particles in the 
fractal which form the defects at the fractal periphery. The evolution of the shape of the 
large droplets slows down with the growth of their size due to the decrease of the droplets 
mobility (see Figs. [3d-c). At the stage when only a few large-size droplets remain their 
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FIG. 8. Similar to Fig. [71 size distributions of islands calculated at different stages of the fractal 
fragmentation (see Fig. [3]) for Ef, = 4 kT, Ae = 0.4 kT and Afi = 2 kT. The corresponding 
simulation time is given in the insets to the plots. 



dynamics is governed to large extend by the interchange of peripheral particles from these 
droplets (see Figs. Eli-f). The large droplets diffuse slowly over a surface and may eventually 
merge. The characteristic time scale for diffusion of an entire large droplet is significantly 
larger than the characteristic diffusion time of single constituent particle, and therefore prac- 
tically can not be resolved within the simulation time limit. However, note that this motion 
can also be simulated with a larger time step. The appropriate value of the time step can 
be estimated using Eq. (|8]). 

Figure El shows the slower evolution of the island size distribution as compared to Fig. [71 
Slowing down of the process is caused by the increase of the binding energy Ei, between 
particles within the fractal. Figure [H shows that as in Fig. [71 a Gaussian-like distribution 
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of a large number of droplets arises immediately after the fractal fast fragmentation. The 
maximum of this distribution slowly drifts towards the larger droplet sizes as the smaller 
islands nucleate (see Fig. [8b). Remarkably, that at the later stages (t = 4 s, Fig. [Hli) several 
maxima arise in the distribution. It is worth noting that this feature of the droplet size 



distribution was also observed in experiment 



Another useful quantity for the characterization of surface structures is the ratio between 



the area and the perimeter of the structure {S/ P ratio) [ij]. This ratio characterizes the 



island topology. Thus, the S/ P ratio for a linear chain of spherical particles is equal to 

where do is the diameter of a particle. Note that the S/P ratio for a linear chain is always 
a constant. The S/P ratio for a compact droplet of the radius, Rd, is equal to 

- = ^. (25) 

P 2 ^ ^ 

It can be easily expressed via the number of particles N in the droplet: 

I = ^VJV. (26) 



In this case the S/P ratio increases as vA^ with the growth of the system size. The S/P 
ratio for a fractal consisting of A^ particles should be larger than in Eq. and smaller than 
in Eq. ( 126|) . Let us now analyze the time evolution of the average (S/P) ratio of the system 
during the fractal fragmentation. The (S/P) ratio for a system of A^ islands is defined as 

iS/P) = ]f E (27) 

^^/r i=l Pi 

where Si and Pi is the area and the perimeter of i-th island, and Nf^- is the number of islands 
in the system. The (S/P) ratio is a useful characteristic for the structures morphology, often 



used in experiment 



The dependence of the (S/P) ratio on time calculated for different sets of the model 
parameters is shown in Fig. and Fig. Hf. Curve 1 in Fig. ^ shows time evolution of the 
(S/P) ratio during the fractal relaxation in the case of the relatively small binding energy 
between the particles being equal to 1 kT. The (S/P) ratio in this case rapidly decreases 
until it reaches the minimum value 0.78 nm, i.e. the S/P ratio which is slightly smaller 
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FIG. 9. Time evolution of the number of fragments/nucleation islands on a surface, -/V/r, during the 
fractal fragmentation process (line 1) and during the nucleation process of randomly distributed 
particles (line 2). The calculations have been performed for a 650 x 750 nm^ substrate with 
periodic boundary conditions. Plots (a) and (b) have been calculated at different values of the 
model parameters: (a) Ei, = 1 kT, Ae = 0.2 kT, = 2 kT; (b) = 4. kT, Ae = 0.4 kT, 
Afi = 2 kT. The insets show the morphology of the system at the end of the simulation. 

than the value for a dimer of particles with do = 2.5 nm. Figure H] shows that the (S/P) 
dependencies to large extend follow the dependencies calculated for (Rmax)- 

The performed analysis provides a lot of useful information on the dynamical evolution 
of the system during fragmentation. However, its direct comparison with experimental 
measurements is rather difficult because the calculated distributions vary with time, but 
the experimental measurements are usually performed for stationary (or quasi- stationary) 
systems. Nevertheless the comparison with experiment is possible if the average life-time T; 
of the studied configuration is greater than the characteristic measurement time T^: 

Ti > Tm. (28) 

Here, T; is defined as the characteristic time-period at which an observable characteristic, 
e.g., the number of fragments in the system, changes within the statistical uncertainty, and 
Tjn is the minimal time-period required to perform an experimental measurement. 
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An important characteristic of the system's stabihty, is the total number of fragments 
Nfj. in the system. At the equihbrium Nf^ fluctuates around the average constant value. 
Note, that Nfr may have similar behavior in a so-called kinetically trapped state, or a 
quasi-equilibrium state which is separated from the equilibrium state by an energy barrier. 
The energy barrier between the kinetically trapped state and the equilibrium state may 
be significantly larger than the thermal vibration energy, therefore the trapped system may 
spend a noticeable lifetime in the kinetically trapped state. This life-time can be sufficient for 
experimental measurements and for holding Eq. (l28l) . This means that the quasi-equilibrium 
value of Nfr may come out different for different initial distributions of particles on a surface, 
demonstrating that different evolution paths may lead the system to different final quasi- 
equilibrium states. Below we analyze two examples supporting this hypothesis. 

Figure [9] depicts the time evolution of the number of fragments/nucleation islands, Nf^, 
calculated (line 1) for the fractal having the initial shape as plotted in Fig. [2|i, and (line 
2) during the nucleation process of randomly distributed particles. The total number of 
constituent particles in both cases is equal to 5182. The size of the substrate used in the 
simulation is identical in both cases, equal to 650 x 750 nm^. Figure [9] shows that the inter- 
particle interaction infiuences significantly the system dynamics. Thus, in the case of the 
weak bonding between particles (i.e. = 1 kT, Ae = 0.2 kT), see Fig. [9^ line 1, the fractal 
fragments into ~2320 islands, i.e. most of the particles in the system are bound in a form of 
dimers. Remarkably, that at these model parameters particles randomly distributed over a 
surface nucleate to approximately the same quasi-equilibrium value Nfr (line 2 in Fig. [9|l). 
The insets in Fig. [9^ illustrate the distribution of particles at the instant t = 4 s in the case 
of nucleation and at t = 2.8 s for the fragmentation. Figure [9^ shows that the system can 
evolve from the very different initial states to the same final state. 

The fragments number evolution with time depend on the inter-particle interaction as 
seen from Fig. [9]3, obtained at larger Eh, Ef, = 4 kT, Ae = 0.4 kT. The quasi-equilibrium 
value of Nfr in this case depends on the initial distribution of particles on a surface. The 
inset to Fig. [9)d shows that both systems have evolved in a group of droplets, whereby the 
size of the droplets created from the initial fractal distribution of particles is larger than the 
size of the droplets created via the nucleation. 

Figure [9] shows that for the chosen model parameters the number of fragments in the 
system becomes constant or changes slowly with time at sufficiently large t value. The 
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FIG. 10. Time evolution of the island size distributions calculated for the nucleation process of 
randomly distributed particles (plot a) and for the fractal fragmentation process (plot b). The 
distributions have been calculated for the same values of the model parameters as in Fig. O). The 
initial fractal shape has been chosen the same as in Fig. [2^. 

resulting static or quasi-static distributions of particles can be compared with experimental 
observations. In the cases when the initial distribution of particles on a surface influences 
the final morphology of the system means the system occupies one of the kinetically trapped 
state. Although the quasi-equilibrium kinetically trapped states do not have the lowest free 
energy, they may live for sufficiently long time to perform experimental measurements of 
the system characteristics. The asymptotic behavior of the fragments distribution with time 
is well seen in Fig. [TOl Figure [10] shows the time evolution of the island size distributions 
calculated for the processes depicted in Fig. [9]d. The island size distribution characterizing 
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FIG. 11. S/P ratio distributions calculated after the fractal fragmentation with different sets of 
the model parameters and the corresponding distributions of island sizes. Distributions (a) and 
(6) are calculated with Eb = 3 kT, Ae = 0.6 Eh, Afj, = 10 kT; (c) and (d) with E^ = 4 kT, 
Ae = 0.4 Eb, Afj, = 2 kT. Insets show the results of experimental measurements for silver fractal 
fragments created via annealing (a) and (6), and by adding of oxide impurities to silver clusters 
(c) and (d) Q- 

the period 0-2 s experiences significant variation, while the distribution during 2-4 s is almost 
static, with only a minor change. 

Figure [TT] shows the island size distributions and the corresponding S/P ratio distribu- 
tions calculated for the fractal fragmentation on the 650 x 750 nm^ with periodic boundary 
conditions. The distributions plotted in Figs. ITTk-b have been obtained with the model 
parameters Eb = 3 kT, Ae = 0.6 kT, A/i = 10 kT at t = 4 s, i.e. well after the fractal frag- 
mentation when the system evolves in the almost stationary equilibrium or quasi-equilibrium 
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state. In this case diffusion of particles along the fractal periphery is the dominating process. 
The increased rate of particle peripheral diffusion leads to the faster island rearrangement, 
and the formation of islands of different size, as seen in Fig. [TTb . The insets in Fig. [TTk-b 
show the results of experimental measurements obtained for silver fractals fragmentation via 
annealing at 600 K. The experimentally measured distribution of the silver cluster island 
sizes is rather broad, with the most probable radius of silver islands ~ 25 nm. A close 
value of 23 nm follows from the theoretical analysis. The discrepancy may arise due to the 
thinner branches of the fractal used in the simulations as compared to the ones analyzed in 
experiment. 

Figures [TTb-d show the island size distribution and the corresponding S/P ratios distri- 
butions calculated with = A kT, Ae = 0.4 kT, Afi = 2 kT. The results of numerical 
calculation are compared with the experimental data shown in the insets to Figs. [TTb-d on 
silver fractals grown with the oxidized silver nanoparticles [l^. In the experiment the most 
abundant radius of the silver cluster islands is 18 nm, being in good agreement with the 
results of our calculations as seen from Fig. [TTH . 

Note that the width and the position of the maximum in the calculated distributions 
shown in Fig.[TT]are rather close to the experimentally observed ones while the absolute value 
of the experimental and theoretical distributions differ quite significantly. This happens 
because we analyze the dynamics of a single fractal, while the experimental measurements 
deal with many fractals on a surface. 



IV. CONCLUSION 

We have performed theoretical analysis of the post-growth processes occurring in a 
nanofractal on a surface using the method, which models the internal dynamics of par- 
ticles in a fractal and accounts for their diffusion and detachment. We have demonstrated 
that these kinetic processes are responsible for the shape of the islands created on a surface 
after the post-growth relaxation. 

The suggested theory is general and can be used in studies of the formation and relax- 
ation processes of different nanostructures deposited on a surface. The developed model 
includes three parameters, which are determined by interatomic interactions in the system 
and could in principle be theoretically calculated for each particular case on the basis of the 
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full atomistic approach for the dynamics of a single particle on a surface. The model param- 
eters can also be obtained from experiment and are specific for different types of substrates 
and deposited materials. In the present paper we analyze the fractal dynamics on a surface 
at various values of the model parameters within a wide range of values and reveal the main 
fragmentation scenarios of the system. 

The paper presents a significant advance in the understanding of paths of the fragmenta- 
tion of deposited nanosystems. It opens a broad spectrum of questions for further investiga- 
tions. Thus, it is interesting to explore the link of the model parameters with the structural 
properties (both electronic and geometrical) of the deposited particles and substrates as well 
as their thermal, mechanical, electromagnetic, etc. properties. Thus, for instance introduced 
model parameters can be determined from the molecular dynamics simulations of different 
diffusion processes occurring on a surface. 

In the performed analysis the deposited particles are assumed to be stiff, i.e. without 
any internal degrees of freedom. However, the particle diffusion over a surface may change 
quite significantly when the particle experiences deformation or changes its phase state. 
Accounting for the detail internal structure and dynamics of particles in the context of 
their diffusion is one of the next obvious steps towards the better understanding of the very 
complex process discussed in this work. 

In the present paper we have studied particle dynamics in 2D. Another obvious important 
extension of the model is to investigate the role of the third dimension in the process of 



fractal formation and fragmentation. This is especially interesting 



many examples of three dimensional fractal systems in biology 85 



to do, because there are 
where the dendritic 



shapes are rather common. Understanding of the growth evolution and fragmentation of 
such systems is very important and may have applications in medicine 
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